Extending Quantum Coherence in Diamond 
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We experimentally demonstrate over two orders of magnitude increase in the coherence time 
of nitrogen vacancy centres in diamond by implementing decoupling techniques. We show that 
equal pulse spacing decoupling performs just as well as non-periodic Uhrig decoupling and has the 
additional benefit that it allows us to take advantage of "revivals" in the echo (due to the coherent 
nature of the bath) to explore the longest coherence times. At short times, we can extend the 
coherence of particular quantum states out from T 2 * = 2.7ps out to an effective T2 > 340/is. For 
preserving arbitrary states we show the experimental importance of using pulse sequences, that 
through judicious choice of the phase of the pulses, compensate the imperfections of individual 
pulses for all input states. At longer times we use these compensated sequences to enhance the echo 
revivals and show a coherence time of over 1.6ms in ultra-pure natural abundance 13 C diamond. 

PACS numbers: 03.67.Lx, 76.70.Hb, 76.90.+d, 33.25.+k 



The nitrogen- vacancy (NV - ) centre in diamond is a 
model quantum system with coherence times of millisec- 
onds, nanosecond gate times, and an optical handle to 
allow initialization and readout of single centres [T] . The 
extraordinary coherence time is essential to proposals 
for using NV centres in quantum information process- 
ing (QIP) [2] or magnetometry [21 [4]. However, the long 
coherence time of the NV defect is not immediately ex- 
ploitable; it is necessary to decouple the electron spin 
from unwanted interactions with its spin-based environ- 
ment, that would otherwise lead to a few /is decay time. 
Decoupling techniques - applying periodic control pulses 
to suppress the interactions - provide a well understood 
solution leveraging decades of use in magnetic resonance. 

The dominant dephasing mechanism of the NV centre 
in high-purity diamond is the surrounding spin bath of 
13 C nuclei in the crystal lattice [5]. The NV is an effective 
spin 1 but a static magnetic field lifts the m s = ±1 degen- 
eracy and we can concentrate on the effective two-level 
qubit system of the m s =0/1 states. With the magnetic 
field along the NV symmetry axis, a secular approxima- 
tion for the effective two-level NV Hamiltonian is given 
by: 

H NV = lo s S z + "jP, +SzJ2 Aj ■ & + n dip. , (1) 

i i 

where S/P are electron/nuclear spin operators; Aj the 
hyperfine coupling to the j'th nuclei; and rldip. the dipo- 
lar coupling between nuclei. An electron spin superposi- 
tion state is dephased by time- variations in the S z oper- 
ator from both fluctuations in the external field and en- 
tanglement with the nuclear bath. At room temperature 
there is an additional incoherent decay when we average 
over all nuclear configurations of the initial maximally 
mixed state during the ~ 10 6 experiment repetitions. 



Because these fluctuations are relatively slow, a spin- 
echo (a 7r-pulse with equal delays, r, before and after) can 
reverse some of the evolution. However, a key part of the 
dynamics is that the quantization axis of a 13 C spin de- 
pends on the state of the electron, due to the anisotropic 
form of the hyperfine interaction [6] ; this gives collapses 
and revivals in the electron coherence [7]. Because there 
are no S x and S y terms in the secular Hamiltonian, for 
a free-evolution time r, the unitary propagator for the 
electron-bath system is U e $ — |0)(0| e ® Ub : o + |l)(l| e t^b,i, 
with the nuclear bath propagators Uq and U\ in the re- 
spective and 1 electron manifolds. This can be un- 
derstood as the nuclear bath measuring an electron su- 
perposition |+) e = i (|0) + |1)); the coherence of the 
electron is determined by the distinguishability of the 
nuclear evolutions in the different halves of the super- 
position. It is straightforward to show [5] that for the 
spin-echo sequence with p nuc . — 1, the expectation value 
of the superposition state cr + = |+)(+| decays as (a + ) = 
\ + \Re {Tr ((tWi) f (UiU ))} . If the hyperfine quanti- 
zation axes are the same, then Uq and U\ commute, and 
there is an echo at all times. However, when the quan- 
tization axes are not co-linear, then there are additional 
echo modulations. The initial echo decays on a timescale 

of a few (j,s « Aj. We can suppress this decay 

by rapidly switching between Uq and U\ (with electron 7r 
pulses) so that they effectively commute - as in a Trotter 
expansion. At longer times, there are unique circum- 
stances where there is a full echo revival. In the |0)(0| e 
subspace, where there is no hyperfine interaction, Uq is 
dominated by the Zeeman term; hence, it is the same for 
all 13 C nuclei and Uq — exp(— it(%2 • + Hdip.) ~ 1 
at t = nQ,L, integer multiples of the Larmor period, Ql- 
Since Uq will then factor out of the trace expression, there 
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are full echo revivals limited by Tidip. and any off-axis 
static magnetic field [Tl [8] (Figure [4] (a)). We can sup- 
press this decay too by adding more pulses; however, the 
minimum pulse spacing must be a full Larmor period and 
the pulse spacings must be commensurate to achieve an 
echo. 
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FIG. 1: NV in a diamond lattice and electron ground state 
level diagram. The electron spin (green) is coupled to the 
surrounding nuclear spins: 14/15 iV (blue) and 13 C (red) in the 
lattice. An example with eight pulses of the three decoupling 
sequences used: CPMG; XY-8 and UDD. CPMG and XY 
have the same pulse spacings but the phases of the XY-8 
pulses are alternated. The timing of the centre of the k th 
UDD pulse with N total pulses is given by sin 2 {-Kk/{2N + 2)). 

Multiple periodic control pulses have been used to 
suppress fluctuating fields since the beginnings of mag- 
netic resonance [5]. Its experimental usefulness has been 
demonstrated in the context of QIP for nuclear bath hy- 
perfine noise [THl [UJ and in ion traps . Of late, there 
has been interest in potential improvements offered by 
non-periodic spacing of the decoupling pulses. Uhrig dy- 
namical decoupling (UDD) [f 3 , a variant of gradient 
moment nulling in NMR[Pl], is provably optimal [15] 
in certain circumstances. Specifically, UDD is optimal 
in a time-expansion of the unitary propagator when the 
noise has a sharp high-frequency cutoff. When the cutoff 
is relatively softer, conventional CP spacing is preferred 
[TBI ITT] . Intriguingly, while suppressing unwanted fluc- 
tuations, both approaches can be used to improve the 
sensitivity of NV magnetometers [H . 

The noise model for a central spin with varying 
anisotropic hyperfine couplings is not conventionally 
studied theoretically, although its pernicious effects on 
control fidelity have been noted [IS]. We show experi- 
mentally that for this model, equally spaced CP style se- 
quences perform no worse than the optimized pulse spac- 
ings of UDD. These experimental tests provide relevant 
results for the realistic case where there are losses due to 
pulse imperfections. However, this loss of fidelity from 
pulse errors is dramatically state dependent. For exam- 
ple, states along the rotation's axis are relatively insensi- 
tive. However, for QIP applications the input state may 
be unknown; this necessitates compensated sequences ro- 
bust to pulse errors. We investigated three decoupling se- 
quences (Figure [I]): conventional CPMG, which for even 



numbers of pulses compensates for pulse errors along the 
direction of the rotation axis; the XY family which are 
compensated for all three axes of the Bloch sphere [20] : 
and finally the UDD which is uncompensated for pulse 
errors. 

We used a standard purpose-built confocal microscope 
setup for optical initialization and readout of single NV 
centres in single crystal diamond. We controlled the ef- 
fective two level system (m s = 0, 1) with microwaves 
modulated by an arbitrary waveform generator to pro- 
vide arbitrary microwave amplitude and phase [28] . 



CPMG Sequence 



64 Pulses 
32 Pulses 
16 Pulses 
8 Pulses 
4 Pulses 
2 Pulses 




10 100 
Total Evolution Time (us) 



UDD Sequence 




Total Evolution Time (us) 
XY Sequence 




♦ 64 Pulses 
■ 32 Pulses 
16 Pulses 
V 8 Pulses 
A 4 Pulses 



10 100 
Total Evolution Time (us) 



FIG. 2: Short-time coherence decay for states along the 
pulse rotation axis with (a) CPMG, (b) UDD and (c) XY 
sequences. Fits are to the phenomenological form s(t) = 
0.5 + Aexp( — {^r) k ) where 3 < k < 6. Error bars are propa- 
gated from photon counting statistics. 

When the input coherence is parallel to the pulse ro- 
tation axis we can explore the extension in T2 with the 
number of pulses and the effect of variable pulse spacing. 
However, the UDD and XY sequences still suffer from 
pulse sequence imperfections and there is evidence of an 
additional early decay with more than 16 pulses. The 
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experimental results are shown for the initial short-time 
coherence decay in Figure [2] and the extracted effective 
T 2 's in Table |U 



TABLE I: Effective T2 (/is) extracted from the fits to the 
curves in Figure [2] For 2 pulses CPMG and UDD are equiv- 
alent and there is no identity XY sequence possible. Exper- 
iments with fixed pulse spacing and going out to 1000's of 
pulses indicate 



Pulses 


CPMG 


UDD 


XY 


2 


11.8 ±0.4 


11.8 ±0.4 




4 


26 ±2 


22 ±2 


27 ±2 


8 


51 ±3 


46 ±3 


50 ±3 


16 


100 ±5 


99 ±8 


91 ±6 


32 


190 ± 13 


185 ± 17 


168 ± 12 


64 


340 ± 25 


293 ± 25 


239 ±21 



We posit that the observed performance of the se- 
quences is due to the nature of the non-Markovian spin 
bath inherent to NV centres in diamond. For the NV 
defect at low fields, the aniostropic hyperfine, and asso- 
ciated coherent dynamics from the dilute spin bath is an 
important non-Markovian noise source not easily treat- 
able as a semi-classical fluctuating magnetic field at the 
electron spin [2T]. We validate our assumption via sim- 
ulations of the NV centre and nuclear spin bath, made 
tractable through the disjoint cluster method [5]. These 
simulations provide insight and confirmation of the valid- 
ity of our experimental results [25] . They also confirmed 
the surprising almost linear increase in T2 as opposed to 
the expected ns dependence for a semi-classical field [2T]. 
They also confirmed that the UDD undcrperformance is 
not caused by imprecision in the pulse spacings due to 
the finite AWG clock frequency, where experimentally 
the timings were rounded to the nearest 2 ns |22j . 

An additional complexity in the NV system is the pres- 
ence of the 14 N hyperfine coupling of 2.1MHz. This 
leads to substantial coherent off-resonance errors in the 7T 
pulses. Our single 32 ns, nominally square, pulses have a 
simulated fidelity of 98.6% (assuming a maximally mixed 
14 N state). Any experimental implementation will fall 
short of this due to errors such as fluctuations in the 
microwave power or transient variations in the phase of 
the pulse at the rising and falling edges. For multiple 
pulse sequences where all the pulses are about the same 
axis, these errors add coherently and rapidly destroy the 
effectiveness for particular input states. For our setup, 
beyond 4 pulses, decoupling for states perpendicular to 
the rotation axis was ineffective (Figure [31). However, for 
QIP applications, the input state is unknown and the 
decoupling sequence should be robust to all input states. 

One approach is to make the decoupling sequence itself 
more robust for any initial state: the motivation for the 
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FIG. 3: CP-8 (circles), and UDD-8 (T) with the input state 
perpendicular to the pulse rotation axis. The deleterious ef- 
fect of the coherent pulse errors leads to two time scales for 
the decay: initial coherent oscillations are rapidly damped as 
they are not properly refocussed and then same long time de- 
cay as the parallel state. Pink (T) Simply modulating the 
phases of the pulses as in the XY-8 sequence (+) can over- 
come this. Robust composite pulses (■) provide similar relief 
despite having five times the number of pulses. Also shown is 
the XY-32 sequence (o) demonstrating that we can preserve 
an arbitrary state to much longer times by shortening the 
pulse spacing. 



XY family of sequences [20]. By alternating the phase 
of 180 degree pulses about the ±X and ±Y axes and 
concatenating appropriately the pulse sequence can be 
compensated for pulse errors and the overall errors made 
isotropic (Figure [3). Indeed, our experimental results 
show that this sequence provides useful refocussing for 
both input states where sequences with all the pulses 
about the same axis are not useful. The tradeoff is that 
this sequence will not perform as well as CPMG for a 
known input state (Figure [2]) . 

A second approach is to use composite pulse sequences 
to make the net rotation of individual pulses more robust 
[23]. We chose the 180 30 - 180 - 180 90 - 180 - 180 30 = 
Z 60 1 80 u sequence for its robustness to resonance offsets 
and ease of calibration [3DJ . The additional Z rotation is 
easy to absorb in an abstract reference frame shift. Im- 
plementing this composite sequence with 46ns Gaussian 
pulses (to avoid overlap of transients at the pulse edges) 
gives a ideal fidelity of 99.95%, and as seen in Figure [3j 
an overall improved performance. Combining composite 
pulses with XY sequences provides even better perfor- 
mance. 

The effect of the multiple pulse echoes can be particu- 
larly dramatic when we observe the echo revivals. These 
revivals will be useful in refocussing the electron spin 
coherence while performing nuclear gates in QIP appli- 
cations [2 . The revivals decay due to off-axis fields and 
nuclear dipole-dipole coupling of up to a few kHz. By 
applying more than one echo pulse, then we can sup- 
press the decay of the revivals at the expense of less fre- 
quent revivals: the revivals only occur when the shortest 
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FIG. 4: The effect the robust XY muftiple pufse sequences 
on the revivafs at short times (a) and ionger times (b). At 
short times the periodic revivals of the single spin echo are 
observed to decay with at 220 fj,sT2- With multiple pulses, 
the revival frequency is reduced proportionally but the T2 is 
greatly extended to over 1.6 ms. The envelopes are fit to 
s(t) — 0.5 + Aexp(— (^) 3 ). Since the revivals only occur at 

certain multiples of the bare 13 C frequency, we observe only 
at these times to reduce experimental time. 



pulse spacing corresponds to a nuclear identity operation 
in the m s = manifold. Hence, it is also not possible 
to see echo revivals with the unequal pulse spacing of 
UDD. We were able to demonstrate an over seven-fold 
increase in the "T2" over the single spin echo time by 
using the robust XY-4, XY-8 and XY-16 sequences from 
220/xs to over 1.6 ms (Figure [4]). This newly revealed 
extraordinary coherence time is then comparable to the 
1.8 ms reported for isotopically purified diamond, some 
of the longest room-temperature coherence times for a 
solid state electron spin [2"4"] . 

In summary, we have experimentally demonstrated 
dramatic increases in the effective dephasing time of a 
nitrogen vacancy centre in diamond by using robust se- 
quences and composite pulses to suppress both the fluc- 
tuations that lead to the dephasing and the intrinsic er- 
rors in the pulses themselves. We expect these robust 
sequences will prove useful not only in both NV mag- 
netometry and QIP applications but also in other solid- 
state QIP implementations such as quantum dots [TU] 
or superconducting qubits |25j . The single spin echo co- 
herence time of our diamond was not particularly long 
at 220 fis, as coherence times of up to 600^s/1.8 ms 
in 1.1%/0.3% 13 C have been reported. These sequences 
should similarly enhance the coherence times of these ap- 



proaches and help them approach the ultimate T± limit. 
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SUPPLEMENTARY MATERIALS 
Experimental Methods 

We used a purpose-built confocal microscope setup for 
imaging single NV centres in a commercial (Element 6, 
electronic grade, CVD-grown, < 5ppb native N impu- 
rities) single crystal diamond. The density of centres 
and photon counting rates strongly indicate we were fo- 
cussed on a single NV, athough we did not explicitly con- 
firm this through photon correlation measurements. Mi- 
crowaves from a low phase-noise source (Agilent E8241A) 
are brought to the centre via a 25yLtm diameter copper 
wire mounted on top of the crystal. The microwaves 
were single-side-band modulated at an IQ mixer (Marki 
Microwave IQ0255LMP) by a 100MHz I.F. created by 
a 500MHz AWG (Tektronix AWG5002) to allow arbi- 
trary amplitude and phase control. This I.F. approach 
allows us to leverage the 14-bit vertical precision of the 
AWG to set both the phase and amplitude of the mi- 
crowaves. Non-linearites in the amplitude were corrected 
for by measuring the spin-dynamics at a range of power 
levels and fitting a non-linear correction curve: (er z ) was 
measured as the power was increased from 1 to 100% and 
the Rabi oscillation curve was heuristically fit to create a 
non-linear correction table. An 10W microwave amplifier 
(HD Communications HD20323 ) gave nutation frequen- 
cies of 15MHz. The T| of the centre measured from fit- 
ting a Gaussian decay to Ramsey fringes was 2.7/Lts; the 
Ti measured in a single spin echo was 220/iS and the T\ 
for the m s = electron state was 6 ms. The measured 
fluoresence (averaging 0.015 photons/350 ns counting pe- 
riod for the m s = sublevel) was referenced and scaled. 
Before each individual run of the experiment, references 
were taken for both the bright state (a 3 [is delay af- 
ter optical pumping) and the dark state (using a robust 
broadband adiabatic state-to-state inversion pulse to the 
m s = +1 state). A static magnetic field of a few tens 
of Gauss was applied along the NV axis using a 3-axis 
Helmholtz coil setup. The field was aligned by maximiz- 
ing the Zeeman splitting as measured through CW ESR 
spectra and comparison of the revival frequency in the 
spin-echo experiment with with the Zeeman shift of the 
electron resonance frequency. This static field splits the 



ra s — ±1 levels; we concentrated on the m s = 04+1 
transition at 2.98 GHz. 

Disjoint-Cluster Simulations 

As shown in the main text, the electron spin dephas- 
ing is determined by the distinguishability of the nuclear 
spins' evolution. Simulating this enormous Hilbert space 
of 100's of nuclear spins can be made tractable by break- 
ing the system up into clusters of strongly coupled nu- 
clear spins: the so-called disjoint-cluster method. The 
code was written in Matlab and is available from the au- 
thors. 10 lattice sites in every direction from the defect 
were considered and 13 C spins were randomly distributed 
at 1.1% giving approximately 800 13 C spins. The nu- 
clear clusters were considered with up to 6 13 C spins 
and the full second order perturbation Hamiltonian ac- 
counting for electron mediated nuclear-nuclear interac- 
tions and the hyperfine-enhanced nuclear Zeeman effect 
was calculated. The simulations only consider perfect 
zero-time pulses and so we can only explore the effect of 
pulse spacing: CP equal time spacing versus UDD opti- 
mized spacing. 

Not surprisingly, the specific form of the decay curves 
varied significantly from run to run, depending on the dif- 
ferent random distributions of 13 C spins around the sim- 
ulated defect. The sample curves shown in Figure [5] all 
correspond to a simulation with a particular defect with 
no stronghyperfine interactions (> 300 kHz). This is sim- 
ilar to the experimental situation where certain centres 
are avoided because of their unresolved hyperfine inter- 
actions lead to a rapid decay and broad CW lines. The 
simulations confirmed the experimental CP-style equal 
pulse spacing out-performs or equals the performance of 
UDD, the scaling of T2 with the number of pulses and 
also showed remarkable agreement with the experimen- 
tal results. 

TABLE II: Effective T 2 {/is) extracted from the fits to the 
simulated curves in Figure [5] 



# of Pulses CPMG UDD 

2 11.1 ±0.3 11.1 ±0.3 

4 25 ± 1 20 ± 1 

8 45 ± 1 46 ± 7 

16 103 ± 3 95 ±3 

32 187 ±7 189 ±7 

64 374 ± 16 374 ± 14 



We can also use the simulations to look at the long- 
time echo revivals and how they are affected by multiple 
pulses. Figure [6] shows the revival at 32 Larmor periods, 
the longest we looked at experimentally. It shows both 
the increase in fidelity with increasing numbers of pulses 
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FIG. 5: Sample decay curves for a particular 13 C distribu- 
tion from the disjoint cluster method for (a) CP style pulse 
spacings and (b) UDD pulse spacings. 



as expected, but tends to overestimate the observed co- 
herence times. Perhaps this is due an additional decay 
mechanism though other nearby defects or a slightly mis- 
aligned magnetic held. 
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FIG. 6: Sample revival curves for multiple pulse echos 
showing the expected improvement with additional pulses. 
Whereas the single pulse echo has decayed to almost nothing, 
16 pulses can create an almost full revival. 



rather than fixing the number of pulses and increasing 
the time. This allows us to explore the limit of our de- 
coupling approach for certain input states. The results 
show that at short pulse spacings we get coherent oscilla- 
tions in the decay as we rapidly modulate the dynamics 
but we extend the coherence time out to greater than 800 
/is. At longer pulse spacing we get a more conventional 
exponential decay as there is sufficient time for small but 
irreversible dynamics in the nuclear bath. At even longer 
pulse spacings there is no high-fidelity echo for even one 
pulse due to the rapid ESEEM decay. 





1.1 




1 


as 


CO 


0.9 










m 


0.8 






o 


0.7 


O 




a. 


0.6 




0.5 




200 400 600 800 1000 

Total Evolution Time (us) 




200 400 600 800 1000 

Total Evolution Time (u.s) 




200 400 600 800 1000 

Total Evolution Time (us) 



FIXED PULSE SPACING CPMG 

The conventional approach to CPMG from NMR is to 
fix the pulse spacing and increase the number of pulses 



FIG. 7: CPMG decay with fixed pulse spacings. r refers to 
the initial delay in the (r — n — 2r — n — r) n sequence. Fits 
are to an exponential decay but the quality of the fits is not 
high enough to give a decay constant. 



